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fS| ■ We investigate the possibility that Uncar arrays of atoms can guide matter waves, much as fiber optics guide hght. 
; I ' We model the atomic line as a quasi-lD array of s wave point scattcrers embedded in 2D. Our theoretical study 
^ , reveals how matter wave guiding arises from the interplay of scattering phenomena with bands and conduction along 
' the array. We discuss the conditions under which a straight or curved array of atoms can guide a beam focused at 
one end of the array. 

I. INTRODUCTION 

In this paper, we discuss the possibility of using a line of atoms to guide matter waves, such as electrons or other 
atoms. Periodic arrays often behave as waveguides; perhaps the most familiar example is electrons propagating in a 
metal. Such waveguides can be engineered as well, such as the guiding of electromagnetic waves in a photonic crystal 
comprised of aluminum rods (see e.g. [H). Although the problem of scattering from periodic arrays is an old one, 
it arises now in a completely new context, as recent technology allows structures to be engineered from individual 
atoms. Single chains of Au atoms, for example, have recently been deposited on NiAl(llO) as well as on Si(553) 
Q . Although our formalism does not treat this situation specifically, artificial arrays of atoms also arise when atoms 
are confined to individual sites in an optical lattices. One might also imagine the individual scatterers being an array 
of quantum dots, which could be used to guide electrons. 

We model the atomic array as a set of s wave scatterers. The question of relevance from an experimental point of 
view is whether we can use such an atomic array to guide waves. Theoretically, this question translates into whether 
eigenstates of the system exist which are propagating along but are evanescent transversely to the array. In order 
to search for these states, we employ scattering theory. Along the array, the system's periodicity gives rise to Bloch 
waves, band structure, diffraction, and other features familiar from solid state systems. However, as the array is finite 
in the transverse direction, it acts like a partially transparent wall in an unusual interplay between Bloch waves and 
scattering theory: transmission and reflection coefficients, along with resonant phenomena for propagation through 
the wall coexist with conduction along the quasi-lD array. Although the same coexistence exists for any real finite 
sample of periodic material, it is especially evident and exposed in the quasi-lD periodic array. 

In order for quasi-lD guiding to exist, the individual scatterers making up the wall must be attractive. However, 
assuming elastic scattering, a freely propagating incident mode can be captured and guided along the array only if 
the array supports conducting states a,t E > 0. It is easy to motivate the fact that that an array of atoms (with each 
atom supporting a.n E < bound state) can give rise to an overall E > conducting state. 

Consider a double well potential, such that each individual well, when taken alone, admits a single bound state. At 
large separation, the wells will give rise to a degenerate doublet of symmetric and antisymmetric bound states. As 
the wells come closer together, the antisymmetric state rises in energy until it is pushed to i? > and becomes a p 
wave resonance, as described in Ref. @, which is closely related to the concept of proximity resonances Q. In Fig. 
[1] the bound state energies for two and four scatterers are shown; as the scatterers come close together the highest 
energy bound state is pushed above threshold. For the infinite wall, one way to understand the presence of i? > 
guided states for the wall is to consider the limit where the scatterers are infinitely close together, thereby effectively 
forming a potential trough of some width w: 

^(^^y)-]^ \x\>w. 

Since the above potential is separable, if the effective ID well (about the x direction) admits an E < bound state, 
the trough admits a continuum of E' > states which are bound along the x direction, because any wavenumber can 
be taken in the free y direction. Thus by analogy to other, more familiar systems, it seems reasonable that a wall of 
attractive scatterers could behave as a waveguide. 
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FIG. 1: Bound state energies for two scatterers (left) and four scatterers (right) separated by distances d. The dashed line is 
the bound state energy of an individual scatterer. For large d, the ground state is the completely symmetric eigenstate, while 
the highest energy state is the completely antisymmetric eigenstate. As the scatterers come closer together, the antisymmetric 
eigenstate is pushed to positive energy. The scatterers are point scatterers modelling the s wave scattering from a cylindrical 
well of depth Vb — —0.8 and radius a — 0.2. 



We further note that at low energies, a curved chain of discrete, attractive atoms approximates a continuous, curved 
waveguide, suggesting an efScient numerical method for modeling continuous waveguides via multiple scattering 
methods. A related method is the boundary wall method in 

The remainder of this paper is divided into five parts. We begin, in Section |TT1 with a brief review of Foldy's 
method of multiple scattering. In Section Hill we apply Foldy's method to solve the Lippmann-Schwinger equation for 
a general periodic array of scatterers. In Section HVl we investigate diffraction and threshold resonances, uncovering a 
type of quasibound states which are related to guiding. In Section|V]we search for true guided states, and we calculate 
the band structure for a single infinite wall of attractive scatterers. We demonstrate that when individual scatterers 
have E < bound states, an array can support conducting states at positive energies, E > 0. Finally, in Section IVTl 
we numerically demonstrate that if the array's symmetry is somehow broken, these conducting states can be used to 
guide a beam; i.e., the array can capture an incident wave focused on one end of the array, forcing it to propagate 
along the array and emerge on the other end. 

II. BACKGROUND: FOLDY'S METHOD 

The physical underpinning of all the resonance and interference phenomena which wc discuss in this paper is 
multiple scattering; both within each unit cell and between different unit cells. To model multiple scattering, we make 
extensive use of Foldy's method ^S], which we briefly review here. Consider a wave 0(r) incident on a collection of 
N identical point scatterers at positions {fi, . . . ,r;v}, where fn = (xmUn)- Applying the t matrix formalism we 
characterize a single scatterer at by its t matrix, 

t = s{k)\r,){r,\ 

where s{k) is a function of the wavenumber, k = \J1mE jh. The functional form of s(fc) is chosen to simulate the 
scatterer of interest under the constraint that s(fc) must satisfy the optical theorem 

ofc 2 

Ims(fc) = |s(fc)|'. (1) 

m 

In the t matrix formalism, the Lippmann-Schwinger equation for multiple scattering is 

N 

^{f) = d^ir) + s{k) Mn)Go{r, (2) 



where 

9771 r i 

(3) 



--H','>{k\r-f,\) 
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is the 2D free-space retarded Green's function satisfying 



(V2+fc2)Go(f,fo) = ^<5(2)(f,fo), 



and the various ipi{fi) in Eq. ([2]) are defined recursively as 



N 

Mn) = (t>{n) + s{k)J2^J3{rj)Go{n,rj). (4) 

tpi{ri) is the effective incoming wavefunction evaluated at the i^^ scatterer after scattering from each of the other 
scatterers, excluding the singular self-interaction of the i*'' scatterer. In the remainder of this paper, we set h = m = 1 
and denote s{k) as s. 

Defining two x 1 column vectors, (j) and ip , whose i^^ elements are given by = (j){ri), ipi = V'('^i): the Lippmann- 
Schwinger equation can be written as a simple matrix equation by inverting Eqs. ([2])- ([11) to yield 



where 



and the matrix G is defined by 



tf=M-^$ (5) 



M = 1 - sG 



G excludes the singular self-interactions of each scatterer with itself. 

Substituting the values of ipi from Eq. ([5]) into Eq. ([2]) yields an expression for the scattered wavefunction: 



N 

^{f) = 0(f) + s ^ Go{r, n) (m- V) ^ • (7) 



III. MULTIPLE SCATTERING FROM GENERAL PERIODIC STRUCTURES: CONNECTION TO 

BLOCK WAVES 

The problem of guiding is related to finding conducting eigenstates for an array of scatterers. We discussed in 
Section U how Bloch waves and scattering phenomena coexist in such systems. The usual 3D approach of reciprocal 
lattices could be adapted to 2D. We instead apply multiple scattering theory, considering all of the multiple scattering 
events within each unit cell as well as between different unit cells. The multiple scattering approach, although more 
involved, yields detailed information about interference processes and also generalizes more easily to the introduction 
of disorder into the lattice. 

In this section, we apply Foldy's method to solve the Lippmann-Schwinger equation for a plane wave scattering 
from a general, infinite periodic array of clusters of point scatterers with a Bravais lattice spanned by dy (Fig. [2]). 
The multiple scattering solution in Eq. ([5]) would seem to require inversion of a biinfinite matrix. In this section, 
however, we reduce the solution to inversion of an iV x iV matrix, where N is the number of scatterers per unit cell, 
and demonstrate how Bloch waves arise naturally from multiple scattering theory. The resulting Lippmann-Schwinger 
equation resembles Eqs. (O-® but with an effective scattering strength s and an effective Green's matrix G which 
account for multiple scattering between unit cells. 

Our approach turns out to be related to the Korringa-Kohn-Rostoker (KKR) method [l^, [ll| , with the variation 
that we have represented our scatterers by t matrices and have begun with an s wave approximation. A related 
approach has been applied to study scattering from two-dimensional periodic slabs of scatterers embedded in three 
dimensions [T^ . 
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FIG. 2: Plane wave incident on array with Bravais vector dy. The unit cell is indexed by q, and each individual scatterer is 
indexed by n. 



A. Multiple Scattering from a Periodic Grating 

Denote by fi?'' the position of the n"^ of N scatterers in unit cell q. Further suppose that all N scatterers in a 
unit cell are identical with t matrix given by s (the arguments in this section can easily be generalized to the case of 
nonidentical scatterers). Applying Foldy's method, we can write the Lippmann-Schwinger equation as 

oo N 

^(f) = 0(f) + s J2 J2^o{r,fi'^)^i'\^n'^) (8) 

q— — oo n—1 

oo 

V'i'^(rt^) = 0(^1'^) + ^ E E G„(ri^),rt))^(„^)(r<„P)) (9) 

p= — oo m=l 

where Go{kr) is the free space Green's function {h ~ m — 1) from Eq. Once again, if the effective incident 

wavefunction amplitudes at the scatterers. ipn^ f^'^'' ) in Eq. ([9|), are known, the full wave function is determined 



using Eq. 

It is straightforward to recursively sum and reindex Eq. ([9]) to show that for an incident wave (normalized to unit 
flux along the x direction) : 

(j>{r) = _^e'*^-^ 



the solutions to Eq. ([9]) are reduced to finding the wave function amplitudes for a single unit cell, since 

^(')(r<'')) = e^^''«'^Vi°^(rl°)). (10) 

Let us focus on the N scatterers in the unit cell indexed by q = 0. Using Eq. (fTU|) on the right hand side of Eq. ([5]), 
evaluated for q = 0, we obtain 

oo Af 

^i°n-t^) = ^(^r) +^ E E Go(rl^),rl"))e^'^--^''^(°)(r<?)) (11) 

p— — OO m — 1 

(m,p)#(r«,0) 
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We wish to separate off the m = n term, which corresponds to multiple scattering between each scatterer and its 
periodic counterparts in other unit cells. It is in this term that we must exclude the self interaction, which corresponds 
to m = n and p = 0. Breaking up the sum, we find 



P— — 00 



N 

m— 1 |_p— — oo 

rn^n 



N 



(12) 



In Eq. p2p . Gr is a scalar quantity independent of the configuration of the unit cell and is given by 

Gr = sY,Goik\p\d)e'''yP'^ 



and G is an X iV matrix defined as 



where we used the lattice sum 



m = n 

G{rir'' - fin) m ^ n 



(13) 



(14) 



(15) 



A rapidly converging expression for G{r) can be found in Eq. of Appendix lAl 

Defining vectors of wavelets 



■(0) 



we can solve for the wavelets in the q — Q unit cell: 

where we have defined 
and 



1 - sGr 



M = I SG. 



(16) 



(17) 



(18) 



1 - sGr 

The remaining wavelets in other unit cells are then simply determined by Eq. (fTO 



(19) 



■(9) _ ^ikyqd^iO}^ 



(20) 
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The full wavcfunction is finally given by substituting Eq. p?)) and Eq. pO|) into Eq. ([8|): 

oo N 

V;(f) = ^{r) + S ^ ^Go(r,riP))e*'^'''P^(M-V^°') (21) 

p^ — OO n— 1 ^ 

N oo 

= (/)(f) + 1^ (m-10(o)) ^ Go(f-7l''\p#)e''=«P'* (22) 

n— 1 ^ p— — oo 

N 

= 0(f) + S^G(r-r<„o))(M-i0^°^) (23) 

Note the similarity between Eq. ([^5]) and Eq. Q: one can go from a single unit cell to a repeating array simply by 
replacing the t matrix s{k) with its renormalized version s{k), and the free space Green's function Go(r) with the 
effective Green's function G(f). The wavefunction in Eq. ((2T|) is a Bloch wave since 

V'(r + dy) = e^'^^'^tpif). 

We note that the renormalization and interference effects in a periodic grating are very similar to the effects encoun- 
tered when a scattercr or a cluster of scattcrers is placed in an external confining potential (sec e.g. [a, [11). The 
similarity arises because scattering in a confined geometry is also effectively a multiple scattering problem: A particle 
can scatter once from the target, reflect from the confining potential, and scatter again. In the case of a cluster of 
scattcrers confined to a hard walled or periodic waveguide, the mapping to an array is in fact explicit; applying the 
method of images, the cluster becomes an infinite array where the effective "unit cells" are images of the confined 
cluster of scattcrers. 



IV. DIFFRACTION AND QUASIBOUND STATES 



The result in Eq. ((23)) is in terms of superpositions of spherical waves. Alternatively, Eq. ([23]) could be written 
in a basis of plane and evanescent waves. Substituting Eq. (|A3|) into Eq. ()23p yields a plane plus evanescent wave 
expansion for the scattered wave: 



^/)(f) 



IS 
1 



oo N 



g— — oo 71 — 1 



1 ,fcl«)|.-.W|^,fc(,)(^_^(0)) 



where the wavenumbers of the diffracted beams are quantized by the Bragg condition: 

Uq) _ l(0) , ^ 
Ky - Ky + ^ 



Uq) - 

A/™ - — - 



(24) 

(25) 
(26) 



The values of q for which k!^^ is real correspond to diffracted plane waves, while the remaining values of q correspond to 
evanescent waves. The set Q of open channels, corresponding to diffracted plane waves, is defined by Q = [Qmm, Qmax] 
where 



(fc + ky)d 

(k — ky)d 
"hi 



(27) 



(28) 



Eq. (|24p is in essence Bragg diffraction: the scattered wave consists of a finite number of diffracted plane waves 
propagating at the Bragg angles given by Eqs. (P5)) - (|26p and an infinite number of evanescent waves. 

In the far field, defined by \x\ » |a;„| for all n, only the diffracted beams survive, and the transmitted and reflected 
wavefunctions are 



q&Q 



(29) 
(30) 



qeQ 
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(b) 




2jr 43ir fi>r llh^ Hn- I4jT lf>x 



2^1. A7t 6JT Stt 10 jt I2jr 14 jr 16 jr 



FIG. 3: (a) Transmission coefficient for scattering of an incident beam with wavenumber k from a waU of hard disks with radius 
a. The incident beam is at normal incidence (fcy"' =0). At low energies, the beam is fully reflected; at high energies the beam is 



almost entirely transmitted. Note the transmission resonances at = fc}," . (b) Structure of the resonances at a = 0.1 (dashed 
line in (a)). This figure also illustrates the breakdown of the s wave approximation at high energies; classically, the limit of the 
transmission should not be exactly 1. 



where the reflection and transmission coefficients, i?^ and T^, are given by (using Eq. (|24p ') 



Ra 



(31) 
(32) 



The quantities \Rq\^ and iTql^ correspond to the probabihty of reflection or transmission into the mode; the q = 
mode is the specular component. The total reflection and transmission probabilities for the wall of scatterers are 
given by 



qeQ 



(33) 
(34) 



The value of renormalized t matrix, s, is constrained by combining Eq. p3p and Eq. (j34p with the unitarity 
requirement R + T = I. This constraint is an analog of an optical theorem for the grating. The optical theorem for 
the single wall of scatterers is derived in Appendix [Bj 



A. Threshold Resonances and Quasibound States 



Since the signature of a bound or quasibound state is often a transmission resonance, we begin by examining the 
transmission coefficient of the array, Eq. As illustrated in Fig. [3] for the single chain of atoms, resonances occur 

at energies such that ki^^ —>■ 0. Such resonances, called threshold resonances, correspond to energies at which one 
of the evanescent waves becomes a propagating diffracted beam. These threshold resonances are purely due to the 
periodicity of the array, and the resonance energies are independent of the type or configuration of the scatterers within 
an individual unit cell. Threshold resonances occur in atom-surface scattering as selective adsorption resonances [lH , 
in X-ray diffraction as emergent beam resonances, and in acoustics as Parker resonances [iGj . 
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Near a threshold resonance, i.e., as fc — ^ fey''' ± e, the lattice sum in Eq. (jA7p yields 



lim Gr = -^—^ (35) 



From Eq. (0 and Eq. 

1 



lim S = ^ = id\l±2k^/'e. (36) 



From Eq. p6p . the dependence of s on the bare t matrix s{k) cancels entirely near the threshold resonance. Exactly 
on threshold, s = 0, which from Eqs. (j3ip - ((32|) yields T = 1; the array becomes entirely transparent, and the incident 
beam is entirely transmitted (although with a phase). This transparency, which is a form of the Ramsauer-Townsend 
effect, is consistent with flux conservation, as the guided beam carries no flux away from the array. 

Inserting Eq. ((36)) into Eq. p4)) . the scattered wave near threshold consists entirely of the emergent beam traveling 
along the ij axis: 

N 

hm V(r) = e*^>y fM-V<"^) e-V24"H-4'"le^fc^"te-?^i°>) + 0(ei/2) (37) 

i^ — Ky e n— 1 

N 

lim 



V;(f) = e^^'-'> ^ (m- e^V^|--S^'le''=l'''fe-.<">) + o(ei/2) (38) 



Fig. m illustrates the probability density above (Fig. IDJb)) and below (Fig. IHa)) a threshold resonance for a wall 

of hard disks. For clarity, the incident wavefunction e^^'^ is omitted from 'ip{r) in Eqs. ([37|l -(|38 p when plotting 
|'0(r)P in Fig. U) Approaching threshold from below, the evanescent beam, which is about to emerge, dominates the 
scattering. The scattered state consists almost entirely of a state which is evanescent along the x direction (Fig. 0^) 
but becomes progressively more weakly bound as we approach threshold. At threshold, the scattered state merges with 
the continuum, and for wavenumbers just above threshold, the scattered state is weakly unbound in the x direction 
(Fig. Wp). The threshold resonance thus corresponds to quasibound states which conduct along the wire. Although 
these quasibound states constitute, in some sense, a form of guiding, they are not truly conducting states. The 
0(e^/^) term in Eqs. (|37|) - ([38)) is due to coupling to other unbound states. The 'ip{r) states in Eqs. ((37)) - (|38p are thus 
quasibound rather than being truly bound, and they have finite lifetimes. Similar states, known as Rayleigh-Bloch 
waves, exist in many other physical systems ranging from ocean coastlines [l7j to acoustics [l6j . 



V. CONDUCTING STATES 



Motivated by the presence of quasibound states, we now refine our search to find genuine conducting eigenstates. 
An array of attractive scatterers, in the limit where the scatterers are closely spaced, should resemble a quasi-lD 
potential trough embedded in 2D and give rise to a set of states which are purely bound along the array. Such 
conducting states would be evanescent in the x direction but free in the y direction and thus correspond to states 
with wavevectors 

k = iKxX + kyij (39) 

In contrast to the quasibound states in Section llV Al the truly conducting states would contain only such wavevectors; 
in a time dependent sense; a wavepacket injected into such a state would conduct forever. Unlike the quasibound 
states, bound states depend strongly on the properties of the grating and do not necessarily exist for an arbitrary 
periodic grating. For example, a periodic array of repulsive scatterers, while possessing infinitely many quasibound 
states, would not have any truly conducting states. 

For a single point seatterer in free space, two criteria characterize a bound state: the state is (1) localized (has 
negative energy) and (2) is a pole of the t matrix (the interpretation of this criterion is that the bound state is a 
scattering state which exists even in the absence of an incoming wavefunction). For an array of scatterers, Foldy's 
method transforms the Lippmann-Schwinger equation into the matrix equation. Eq. ([T7]), which can be rewritten as 

(1-sG,)m1"V^°^ = ^^"^^ (40) 



FIG. 4: Threshold resonance for a unit ceU with five randomly placed scatterers (hard disks with a = 0.1) depicted by the dots 
in the figure). The full wavefunction just above and below threshold is given by Eqs. (|37l I38p . For clarity, we plot only the 
scattered wavefunction; the incident wave is a plane wave and would be added to the scattered wavefunction to obtain the full 
wavefunction. (a) The quasibound state just below resonance, at wavenumber kd = 2n ~ 0.3, becomes unbound (b) just above 
resonance, at kd = 2tt + 0.3. 

The existence of a scattered wavefunction for zero incoming wavefunction implies the existence of a homogeneous 
solution to Eq. ([iU)) . States which exist in the absence of an incoming wave must exist at values of k which are roots 
of the secular equation 

(1 - sGr)^dctM = 0. (41) 

This secular equation is strongly dependent on the configuration of the unit cell and must be solved for a particular 
grating. We shall here focus on the single wall of atoms, which is the most common experimental setup QjSB- For 
the single wall, M = I, and Eq. (|4T|) simplifies to 

1 - sGr = 0, (42) 

implying, from Eq. (|19p . that a bound state corresponds to a pole of the renormalized t matrix, s. 

In order to find conducting states, we evaluate S on a grid of values of {kx, ky), and numerically search for poles of 
the form given in Eq. I|39p with corresponding energy 




Bound states with E > correspond to those poles of the renormalized t matrix for which \ky\ > \kx \ ■ We used 
the t matrix in Eq. (jC3[) to simulate an infinite wall of attractive cylindrical wells with depth Vq = —8.0 and radius 
a = 0.2, so that an individual well allowed a single bound state just barely below threshold (at E = —0.05, or 
K = —\/2E ~ 0.32i). Fig. [T] depicts how arrays of two or four such wells support bound states above = as the 
wells are placed close together. Figs. [5ja,c,e) show the manifold of poles in k space for the infinite wall of scatterers 
for three different lattice constants. In each figure, the solid curves correspond to the manifolds of conducting states 
for the array, and the dashed curves correspond to the bound state energies of a single scatterer, given by 
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FIG. 5: (a,c,e) show the curves along which s{iKx, ky) has poles (the symmetry across fc^ = is due to reflection symmetry of 
the array along y = 0). These curves correspond to manifolds of conducting states bound in the x direction, and propagating 
along the array. The dashed curves, for comparison, are contours of wavevectors corresponding to bound state energy of a 
single scatterer. (b,d, f) show the band structure {E vs. ky) for the conducting states. (a,b) correspond to a lattice constant 
d = 3.0; (c, d) to d =1.0, and (e,f) to d = 0.202. As the wells approach each other, the bands are increasingly perturbed from 
the single scatterer limit. In all cases, conducting states are present at iJ > 0. 

A more traditional view of the band structure for each potential is presented in Figs. [UJb, d, f), where the energies 
of the conducting states versus ky have been plotted. Figs. [5ja,b) correspond to a large well separation {d — 3.0), 
Figs. ElJCjd) to an intermediate well separation [d = 1.0), and Figs. [5je,f) to a small well separation [d = 0.202). 
When the wells are far apart, the collective bound state is near the bound state energy of a single scatterer, but as the 
wells approach each other, bound states occurring at _E > start to appear. These states arc particularly interesting 
since they could, in principle, be accessed by a propagating incident wave. As a demonstration. Fig. [5] illustrates 
that an incident propagating plane wave can couple into a state which is bound along the waveguide. Note that the 
periodicity of the guided portion of the wave is not the periodicity of the lattice. 
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FIG. 6; Plane wave at normal incidence (k = 0.7) scattering into a conducting state; scatterers (depicted by dots in the figure) 
are soft disks with Vo = —8.0, a — 0.2 and d = 0.8 (the same scatterers as in Figs. [T]and[5}. This figure is a numerical 
demonstration that a propagating incident wave can couple into an i5 > conducting state of the array, if the translational 
invariance of the array is broken. Coupling into the conducting state with a plane wave does not constitute guiding, as the 
incident wave already has amplitude everywhere. 



We now apply the arguments of the previous sections to demonstrate that a wall of atoms might be used to guide 
waves. By guiding, we mean that if an incident beam is focused on some part of the grating, the scattered wave 
propagates along the grating for a long distance. In Section IIV Al and Section |Vl we discussed the existence of 
states that are either quasibound or bound along the infinite array. These states, as such, are not an example of 
transport; being translationally invariant along the array, they can only be accessed by initial conditions that are 
already translationally invariant, such as plane waves. Transport along the array requires that the array have a 
symmetry breaking point where an incident wave can be injected. In particular, in order for a localized beam aimed 
at the array to conduct, the array must have a defect or an end. The asymmetry will of course affect the properties of 
the grating, e.g. diffraction, impedance, etc. While we can use our results on the infinite array as a basis for studies 
of guiding in related systems, the infinite array itself is not a waveguide. 



where g{9k) — e"'^^''/™)^ . Eq. represents a beam incident from the right, focused on the scatterer at 0, and 
rotated by an angle /3 from the positive x axis. An incident beam [Eq. (|43[) with /3 = (normal incidence)] is shown 
in Fig. [7ja). As w —f 0, (f>{r) approaches a plane wave, and for w — > oo, 4>{f) becomes a spherical wave with its 
form related to the Fourier transform of the spherical wave Jo(fcr). In the following, we are interested in intermediate 
values of w, for which Eq. (^5]) represents a focused beam. 

Using our results for the infinite array as a guide, we take a purely numerical approach to the study of guiding in 
finite and semiinfinite arrays, as these systems are difficult to treat analytically. Analytical studies (e.g. Rcfs. 
[l9| ) typically conclude that scattering from a semiinfinite array reproduces certain major features of an infinite array: 
resonances, diffraction, etc. The main difference is that in the semiinfinite case, a spherical wave emanates from the 



VI. ARRAYS OF ATOMS AS WAVEGUIDES 



A. Guiding by a Semiinfinite Array 



Consider an incident beam 




(43) 
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FIG. 7: Guiding of a focused beam with k — 0.7. The scatterers are indicated by dots. These plots show the probabihty 
densities of (a) the incident beam, focused on the end of the array (b) The scattered wave, which is guided along the array, 
and (c) the full scattered wavefunction (incident plus scattered). The individual scatterers are attractive wells with the same 
form and lattice constant as in Fig. [3e)-(f). 

end of the array. Most studies of semiinfinite arrays begin with the infinite solution and derive this spherical edge 
wave as a correction term. 

We demonstrate guiding numerically in Fig. [71 The scatterers are closely spaced, and the focused beam is aimed 
at the end of the array. In Fig. U^c), a guided state clearly propagates along the array. The guiding occurs at low 
energies, near bound states of the infinite array. Guiding does not occur for repulsive scatterers, confirming that the 
coupling is to the true bound states of the array rather than the quasibound states. 

The guiding is robust and extends to gently curved and finite walls as demonstrated in Fig. [51 where a beam 
focused on one end of a curved wall emerges at the other end. The scatterers are effectively behaving like a light pipe 
for a matter wave. If the wall is curved sharply relative to the wavelength of the incident beam, amplitude leaks out 
as adiabaticity breaks down. 

In this section, we have demonstrated via numerical simulations that the E > conducting states of the array can 
be exploited to guide waves along a semiinfinite wall of attractive, closely spaced scatterers. The physical explanation 
is that a semiinfinite wall of attractive potentials behaves like a trough, and furthermore that because the trough has 
an end, it is possible to couple into conducting states via injecting a beam into the end. 

B. Effective Cross Section for Conduction Along the Wire 

From Figs. [TJSl it appears an incident beam of particles focused on one end of a wall of scatterers can be efficiently 
guided to the other end. In this section, we provide a quantitative estimate of this guiding efficiency. Consider a 
focused beam (f>{f') [Eq. ([43)) ] of energy E = ifc^ with normal incidence (/? = 0) to the wall of scatterers. Such an 
initial incident beam has nonzero fiux ji„c along the x direction, where 

JO 

The scattered flux can be calculated by evaluating the flux of the scattered wave function ipsi^ over any contour 
which completely surrounds the scatterers. Finally, the total cross section can be simply evaluated by integrating the 
scattered flux over the desired contour and dividing by ji„c- For the wall of scatterers, the total cross section, a, can 
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FIG. 8; Guiding of a focused beam along a finite wall; k — Stt. The scatterers (depicted by dots in the figure) are attractive 
disks with Vb = — 10, a = 0.2. These plots show the probability densities of (a) the incident beam, (b) the scattered wave, and 
(c) the full scattered wavefunction (incident plus scattered). In (b), amplitude visibly leaks away near the sharp bend near 
?/ = 0, as a consequence of nonadiabaticity. 



m 



FIG. 9: The "stadium" -like contour used to evaluate the total cross section, a, for the wall of scatterers. The shortest distance 
from any point on the contour to the wall is r„ . For an incident beam focused on the bottom of the wall of scatterers (indicated 
by the arrow), the fraction of particles "transported" along the wall is given roughly by ac^ /f. 



be calculated over the contour shown in Fig. [HI giving 

1 



a = - — / da • lni[-0s (a)VV's(a)] 

Jinc Jc 



(44) 



where (ycuiD) ^^^^ contribution to the total cross section from particles scattered through the upper (lower) semi- 
circles. 

As long the contour encloses all the scatterers, the value of a is independent of the choice of contour. However, 
o'c'u{D) "^'-^^1(2) ^^^^ depend on the value of for the contour shown in Fig. [51 We are interested in estimating 
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FIG. 10: (a) The total cross section, cr, in units of 10*d (wliere d is the lattice spacing) for A'' = 150 square well scatterers 
(Vb = —8, a = 0.2, d = 0.8) as a function of kd. The incident beam is focused upon the bottom scatterer, with u; = 1. Since 
a is much greater than the spread of the incident beam, it can be safely assumed that each incident particles are scattered. 
Therefore in order to obtain the number of particles which are transported along the scatterer wall, one must simply calculate 
acu/f", shown in (b). As can be seen, up to 30% of all scattered particles can be transported along the wire (compared with 
only 0.8% if the scattering was assumed to be isotropic). 

the fraction of scattered particles transported along the wire for an incident beam which is focused on the bottom of 
the wall of scatterers, i.e., we are interested in the ratio of ocujo (assuming each incident particle is scattered). We 
should therefore choose r^, such that r.^j/ L <^ 1, where L = Nd is the length of the wall of N scatterers with spacing 
d. 

Figure [TCTa) gives a (in units of IQ^d) as a function of kd = y/2E for a system of = 150 square well scatterers 
characterized by Vb = —8, d = 0.8, and a = 0.2 (in addition, w = 1 was chosen for cj){r) in Eq. ([15|) ). As can 
be seen from Fig. fTOFa). a increases with decreasing fed, although for the range of kd plotted in Fig. fTOFa). a is 
still much larger than the spatial spread of the incident wave function so that essentially all of the incident particles 
are scattered. In order to estimate the amount of transport along the scatterer wall, the ratio of the scattered flux 
through the other end of the wall is taken with the total scattered flux, which is simply equal to /o'] this ratio is 
shown in Figure [TOT b). Here, the fraction of scattered particles which are transported along the wire can reach up to 
30%, and oscillates with k. Note that the oscillations are related to the actual length of the wall; if the number of 
scatterers (and hence the array's length) is doubled, the period of the oscillations would roughly be doubled too. For 
O.SStt < k < O.SStt, jo oscillates slightly about jo ~ 12%. This means that even for a beam of particles with 
a distribution of incident energies (possibly due to thermal effects), at least 12% of the incident particles will still be 
transported along the wire. Also note that in these calculations, = lOd so that r^jL = 1/40 <C 1. It is interesting 
to observe that for a beam incident only on the lowermost scatterer, and for the case of isotropic scattering (i.e., 
neglecting multiple scattering), the amount of scattered particles transported through the same region relative to the 
total scattered flux would be given by 2ro/(27rL) sa 0.8%. Thus the wall of scatterers can enhance the scattering along 
a given direction by almost a factor of 40. 
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VII. CONCLUSIONS AND FUTURE DIRECTIONS 

In this paper, we have examined scattering from and guiding by quasi-lD periodic gratings of scatterers embedded 
in 2D. This system is a good laboratory for highhghting the coexistence of scattering phenomena, such as transmission, 
reflection, and resonance, with features typical of periodic systems, such as band structure, diffraction and conduction 
along the array. Arranging individual atoms on substrates, often in patterns far more intricate than gratings, has 
become an established experimental technique 0, H, I3j US • The motivation for our study is thus that atomic arrays 
can be built, have been successfully modeled by multiple scattering theories, [111, [l^, [IS] , and can potentially serve as 
a waveguide for other particles. 

In order to investigate the use of atomic arrays for guiding, we have developed a multiple scattering theory (related 
to the KKR method) for quasi-lD gratings of s wave scatterers, embedded in 2D. The central physics is related to 
Bloch's theorem; we can obtain the full scattered wavefunction from the solution for a single unit cell by replacing the 
t matrix of an individual scatterer with its renormalized version and the free space Green's function with an effective 
Green's function. We have used this result to discuss and to examine some general features of scattering from a 
single chain of atoms, such as resonances, quasibound states, and conducting states. Finally, we have demonstrated 
numerically that conducting states of the semiinfinite or finite array can be used to guide waves along straight or 
curved walls with up to 30% guiding efficiency. 

The properties of the scattering and guiding studied here have strong analogs in optical systems, such as planar 
waveguides, optical fibers, and glass tabletops. For the latter, light incident on a sheet of glass at an angle and far 
from any edge is partially reflected and transmitted directly. Some of the light is transmitted and reflected indirectly, 
propagating inside the waveguide with attenuation by incomplete internal reflection. On the other hand, light incident 
on a (symmetry breaking) edge can be partially reflected, but mostly enters the sheet of glass and is trapped inside by 
total internal reflection. Due to its resonant modes, our system has shown analogs to each of these optical phenomena, 
and others. 

This work can be extended in a number of directions. We have chosen to embed our system in 2D because 2D is 
the relevant dimensionality for electrons in surface states scattering from atoms adsorbed on a metallic surface. With 
a change in the free space Green's function, one could, however, revise the entire theory to treat a quasi-lD array 
embedded in 3D, such as might be relevant if the scatterers were atoms confined to an optical lattice. The physics 
of renormalization and interference in a 3D system would be very similar, and presumably identical guiding effects 
would arise. 

Among the reasons that we have taken a multiple scattering approach to the system is that it can be generalized 
to the introduction of impurities or defects into the array. We have discussed how symmetry breaking is required for 
guiding to occur. Introducing an impurity or defect into an infinite array is one way of breaking symmetry, with the 
impurity behaving as an antenna capable of drawing the incident wave into the conducting state. The problem of 
impurities in the array is also interesting in the sense that the array with an impurity becomes a scattering system 
in the conducting mode, yielding, in essence, a scattering theory within a scattering theory. The question of guiding 
in a system with disorder may also be of interest from the perspective of Anderson localization. Interesting questions 
arise of how to maximize the guiding efficiency, for example by impedance matching the array's end to free space. 
Finally, Foldy's method has recently been extended to include Rashba spin-orbit coupling [23|, and it may be possible 
to modify the results of this paper to examine spintronic versions of guiding. 
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APPENDIX A: RELEVANT LATTICE SUMS 



Equation (fT^]) and Equation ([TS]) are in the form of spherical waves. We would hke to find a more useful expression 
for Eq. (|15p . and also a more rapidly convergent expression for Eq. (|13p . Sums involving Hankel functions converge 
to the plane wave limit very slowly. The physical implication of this slow convergence is that edge effects are more 
important than one might expect; many scatterers are required to build an "infinite" wall. 

1. Plane wave form of effective Green's function 

The Green's function for an infinite array of scatterers, as expressed in spherical waves, is 

oo 

G(f)= G'o(f,f„)e''=='"''. (Al) 

n— — OO 

Substituting the two dimensional free spaee Green's function [Eq. ((21)], we find that 

oo 



2 

n— — oo 



^ 

"2^ 



-00 f^'x ^ '- 



Tl — — 00 



where we have used the integral form of the Hankel function in the right half plane (x > 0). Using 



oo 

^ [ 



27r ^ , 2n7r\ 

n— — OO ^ ' 



we can do the integral, yielding the sum of plane waves at real and imaginary Bragg angles, 

oo ^ 

G(r) = -4 V -^e*^'"'"'^^'^"''' 



n— — OO "'3.; 

oo 

e 

n— — oo 



oo ^ 

E -^e-i-l^le-^^ (A3) 



where we have taken the absolute value to ensure convergence, and 

,(„) ^ . 
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define wavevectors oriented at the Bragg angles. Imaginary values of kx correspond to evanescent modes and are 
present only in the near field. 

This idea connects to the phenomenon of a "healing length" when examining reflections from a corrugated wall. 
We note that G{r) is singular as we approach the origin; while the singularity is not evident in Eq. (|A3p , it is readily 
apparent in Eq. (|A2[) . Re-indexing, we find 

The wavefunction is a Bloch wave: 



withi?=ifc2= 



2. Kummer's Method: Extracting the Singularity 

As in Refs. [l3| and [l^l, we want to apply Kummer's method to extract the singularity from Eq. (|A3p and also to 
obtain a rapidly convergent expression for the renormalized scattering strength, s [Eq. (|13p ]. We begin rewriting Eq. 
(ini) as 

Gr = lim [G(f) - Go(f)] 

= \ira[{G{r)-S{r)) + {S{r)-G^{f))] (A5) 
f— >-o 

where S{r) is a sum chosen to cancel the log singularity of the Hankel function. Looking at Eq. (jASp . we choose the 
following form for S{r) : 

f 2mry\ 



n=l ^ ' 



With this choice of S(r)^ we are extracting from the zero energy limit from the singular sum in Eq. (jA4p . evaluated 
(for simplicity) for a normally incident. The choice of Sij) is not unique; our motivation for choosing this particular 
form of S{r) is that we know that it will contain the logarithmic singularity. From Eq. (jAip we see that the zero 
energy limit of the Green's function corresponds to the limit of a Hankel function as we approach the origin, and thus 
the zero energy limit of the sum will be singular. 

We shall proceed to show that although both S{r) and G(r) separately diverge logarithmically at the origin. 



lim[5(r)-Go(f)] 



is finite. We first find a closed form of S(f): 



1 1 

5(r) = -iReV-fe^ 
TT n \ 



a^rl I 2i2ry 



We can now use 



n=l 



to rewrite 
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Using this form, it is simple to take the limit, 



1 / kd\ 1 

hm 5(f) = In — + - In kr. (A6) 



Subtracting the limiting form of the Hankcl function 



i 7 — In 2 1 
lim Go(fcr) = -- + h - In fa 

- - 2 TT TT 



f-+0 



from Eq. (|A6[) . wc find 



limJ5(f)-Go(f)] = -ilnf^^ ■ * ^ 



TT V^T''/ 2 TT 

which is finite. Returning to Eq. (|A5p . we find that another form of G(f) is 



G(r) = -lle'^-l^^l - 1 y f _i_e^^^"'Ne<^-'^+^)^ + J- e^^'kle^^."^)^ - Ae^l-I cos ( ^-^\] 



1 , / kd\ i 7 

In — H ^, 

TT V47ry 2 tt' 

and from this expression, we can calculate Gr- 

Gr = lim [G(f)-Go(f)] 



-i i / 1 ^ 1 { i 7 

kji'd^^ VfcF " 2i |n| TT j " ^ V4^y 2 " ^" 



In the special case where the incident plane wave is normal to the array, this expression is identical to G,. for a 
scatterer at any location in a periodic wire 
now depend on the incoming wavefunction. 



(71) 

scatterer at any location in a periodic wire. In the general case, the value of Gr differs in the values oi kx , which 



APPENDIX B: UNITARITY AND THE SINGLE WALL 

Simplifying Eqs. pT|) - ((32)) to the case of a single chain of atoms, the reflection and transmission probabilities are 

R, = (Bi) 

Tq = 6qQ + Rq. 

The unitarity requirement, T + i? = 1, can be shown to yield the following constraint on s, which resembles the 
ordinary optical theorem, Eq. ([1]): 



ImS = |S|^ y 



Using Eq. (fT5|) and Eq. (jA7p . it is straightforward to verify that s satisfies the unitarity condition, Eq. (|B2p . as long as 
the single-scatterer t matrix s{k) satisfies the free space unitarity condition, Eq. ([1]). In the limit kd — > 0, our results 
approach those for a continuous wall of scatterers. An alternate method for treating a continuous wall is the boundary 
wall method where the wall is discretized into a set of pseudo-scatterers. The boundary wall method has been 
shown to be equivalent to building a wall out of point scatterers-with the exception that the self-interaction, which 
was omitted in Eq. ([5]), has been left in the calculation, rendering the individual pseudo-scatterer unphysical since 
its t matrix no longer satisfies the optical theorem in Eq. ([T|). Our formalism, though significantly more complicated, 
has the advantage of corresponding physically to buiding a wall out of individual atoms. 



19 



APPENDIX C: SIMULATING POTENTIALS WITH A T MATRIX 

1. Hard Disk 

The 2D hard disk is simulated by requiring that the wavefunction go to zero at a radius a from the scatterer. It is 
shown (see e.g. [1^) that this corresponds to a i matrix 

2. Soft Disk 

Another example of a potential for which we can analytically calculate the s wave scattering properties is a soft 
disk of radius a, 

V{r) = Voe{R~a). (CI) 

We are particularly interested in the case Vq < 0, corresponding to an attractive disk. 

The s wave phase shift Sq for the soft disk potential is determined by continuity of the logarithmic derivative of the 
wavefunction at r = a : 

tan(5u - g'^o(^")'^i(g°) ~ kJi{ka)Jo{qa) ^^^^ 
qYo{ka)Ji{qa) — kYi{ka)Jo{qa) 

where 

q = ^2{E - Vo) 



is the wavenumber inside the disk. 

Having calculated the s wave phase shift, it is straightforward to confirm that we can simulate the s wave scattering 
from a soft attractive disk with a t matrix 

m = r^^. (C3) 
1 — i tan do 



The negative energy poles of the t matrix in Eq. (|C3[) occur at 

qKo{Ka)Ji{qa) — KKi{Ka)Jo{qa) (C4) 

where k = \J ~2E. Eq. (jC4[) can easily be shown to be the characteristic equation for bound states of the 2D cylindrical 
well. 



